A study on multi-frequency electromagnetic imaging 



Won-Kwang Park 

Department of Mathematics, Kookmin University, Seoul, 136-702, Korea. 



Abstract 

Recently, a non-iterative multi-frequency imaging algorithm was developed on the based of an asymptotic 
expansion formula for thin, curve-like electromagnetic inclusions and the structure of singular vectors in a 
Multi-Static Response (MSR) matrix, . In this paper, we study the structure of a multi-frequency imaging 
functional and present an improved imaging functional weighted by a given frequency. In order to analyze 
the function structure, we discover the relationship between the imaging functional and Bessel function of 
integer order of the first kind. A set of numerical examples for single and multiple inclusions shows that 
the presented algorithm not only retains the advantages of the traditional function but also improves the 
imaging performance. 
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1. Introduction 

In this paper, we consider the problem of inverse scattering from thin curve-like electromagnetic inclu- 
sion^) embedded in a homogeneous domain. Generally, the main purpose of this problem is to identify 
characteristics of the target of interest, e.g., shape, location, and material properties, but the most important 
information about a crack or thin electromagnetic inclusion problem considered herein is the shape. Conse- 
quently, various shape reconstruction algorithms have been reported in various literature. However, most 
of these algorithms utilize a Newton-type iteration method, which requires additional of an optimized reg- 
ularization term, complex evaluation of Frechet derivatives at each iteration step, and a priori information 
of unknown inclusion(s). However, shape reconstruction via the iteration method with a bad initial guess 
fails if the above conditions are not completely fulfilled. 

For this reason, various non-iterative shape reconstruction algorithms have been developed, such as 
Multiple Signal Classification (MUSIC) algorithm (a, |17J, [20], topological derivative strategy 11X31. 11611 . and 



linear sampling method [7, 12]. Recently, a multi-frequency based algorithm was developed for obtaining 
a more accurate shape of unknown inclusions. Related articles can be found in |1 Eil El EE ES E3 E? 



and references therein. However, studies have applied this algorithm heuristically and therefore certain 
phenomena such as appearance of unexpected ghost replicas cannot be explained. This motivated us to 
explore the structure of multi-frequency imaging algorithm, find optimal conditions for good imaging 
performance, and improve it to obtain better results. 

The analysis of this algorithm is based on the fact that multi-frequency based imaging functional can be 
represented as a definite integration of the square of the Bessel function of order zero and of the first kind. 
Using |[6]l in section[2j we obtain this imaging function at a given frequency to as follows: 



/■ 



Jo {ujx) dw = w I Jo (lux) + Ji (tax) ) + / Ji (lux) dw 
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where x denotes the location of the unknown target. Hence, one of the way of improvement is to eliminate 
the last term of the above identity (although it does not significantly contribute to the imaging). For this 
purpose, we design an alternative imaging functional weighted by applied frequencies. 

The rest of this paper is organized as follows. In section [2] we introduce the two-dimensional direct 
scattering problem and an asymptotic expansion formula for a thin electromagnetic inclusion, and then 
review the multi-frequency imaging function presented in H, [l3,[l^,[l^l . In Section [3) we discover and 
specify the structure and properties of the existing multi-frequency imaging function. We then design an 
improved imaging function weighted by several applied frequencies and analyze its structure in order to 
investigate its properties. In section 2) numerical results are discussed in order to verify the effectiveness 
of the improved imaging function. Finally, section [5] concludes the paper in brief. 

Finally, we would like to mention that although the shape constructed using the proposed algorithm 
does not match the target shape completely, considering it as an initial guess for an iterative algorithm will 
guarantee a successful reconstruction; for details, refer to j2l.l8l.Elll. 



2. Preliminaries 



2.1. Direct scattering problem and asymptotic expansion formula 

We briefly survey two-dimensional electromagnetic scattering from a thin, curve-like inclusion in a ho- 
mogeneous domain. For this purpose, let O and T denote a homogeneous domain with a smooth boundary 
and a thin inclusion, which is characterized in the neighborhood of a simple, smooth curve 7: 

T = {x + pj7(x) : x e 7, P € {-h, h)} , 

where h specifies the thickness of T. Throughout this paper, we denote t(x) and T}{x) as the unit tangential 
and normal, respectively to 7 at x. 

Let £0 and /^o denote the dielectric permittivity and magnetic permeability of O, respectively. Similarly, 
e and /i denot the permittivity and permeability of T, respectively. At a given non-zero frequency uj, 
let (x; uj) be the time-harmonic total electromagnetic field that satisfies the following boundary value 
problem: 

V • ( -x(T) + — x(0\T) J Vu (/) (x; uj) + lo 2 ( e*(T) + e x(«\T) ) u {l) (x; uj) = for x£(! 



Mo 



1 du®(x:w) 1 de luj6 ' x , 

for x e oil 



(1) 



^0 <9n(x) ^0 9n(x) 



with transmission conditions on 9f2 and dT. Here, {61 : I = 1, 2, • • • , Q} is the set of incident directions 
equally distributed on the unit circle S 1 , and \ denotes the characteristic function. Let ujp (x; uj) — e lu} ° 1 x be 
the solution of (Q without T. Then, due to the existence of T, the following asymptotic expansion formula 
holds. This formula will contribute to development of the imaging algorithm. A rigorous derivation of this 
formula can be found in |6|]. 

Theorem 2.1 (Asymptotic expansion formula). For x £ 7 and y e d£l, the following asymptotic expansion 
formula holds: 

u^{y-uj)-4\y-u)^hu^{y-uj) + o(h), 
where the perturbation term (x; uj) is given by 

uW (y; uj) = J [uj 2 (e - e )u { l) (x; w)A(x, y; uj) + Vu« (x; uj) ■ M(x) • VA(x, y; uj)j d 7 (x). (2) 
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Here, A(x, y; w) is the Neumann function for Q that satisfies 
1 



/'-o 



and M(x) z's a 2 x 2 symmetric matrix such that 
• M(x) has eigenvectors r(x) and ^(x), 



-V-VA(x,y;aj) + u%A(xj;w) = -i(x,y) /or xe!) 



1 <9A(x,y;w) 
Ho <9n(x) 



/or x G <9ft, 



• eigenvalues corresponding to t(x) and rj(x) are 2 — awd 2 — respectively. 

2.2. Introduction to multi-frequency imaging function 

Now, we apply the asymptotic formula 10 in order to introduce the multi-frequency imaging function. 
Throughout this section, we assume that for a given frequency oj = 7 is divided into M different 
segments of size of order 4, and only one point, say x m for m = 1, 2, • • • , M, at each segment contributes 
to the image space of the Multi-Static Response (MSR) matrix A(cj) based on the Rayleigh resolution limit; 
for details, refer to [? EH]. 

For introducing multi-frequency imaging algorithm, we consider the decomposition of MSR matrix 
A(lu) — (Aji(uj)) G C Pxl 3, whose element Aji(u>) is the following normalized boundary measurement: 



A,-,(w) 



u\'{y;uj) — ^— — dS(y) = uj 



1 dn dn(y) 
w 2 length (7) 



(e - £ ) + #j • M(x) • 6 t )e-^ ( ^- ei) - x d7(x) 



,2 length(7) '"" 

M 

rn—1 



M 
1 



E 



(e - e ) + 2 ( - - — J ^ • r(x m )0; ■ r(x m ) 
A* Mo/ 



/./, 



Mo M5 



i9, • ri(x m )9i ■ r)(x m ) 



where length(7) denotes the length of 7 and (x; w) = e^ l ^ j x , j = 1, 2, • • • , P, is a test function. Note 
that : Z = 1, 2, • • • , P} is the set of unit vectors equally distributed on the unit circle § . 
Based on the above formula, A(w) can be decomposed as follows: 



A(w) = J(w)B(w)L(cj), 
where B(cj) G jj3A/x3A/ ^ g a diagonal matrix with component 



(3) 



w 2 length( 7 ) ( e-s 



M \ 
and J(w) G C Px3M and L(w) G C 3MxQ are written as 





l(x m ) 



J(w) = ( Ji(w), J 2 (w), • • • ,Jm(w) ) and L(w) = ( Li(w), L 2 (w), • • • ,L M (w) ), 
respectively, where 



J»n(w) 



\ e^'' 1 *, #p • r(x m )e-^- x , P • r,(x m ) 



^•J7(x m )e-^ 1X \ 
*2 ' »7(x m )e-^' 92 - x 



/ 



and 



/ e ^ 0l X -, 0i -T(x m )e^ fll - X , 6»i ■ r;(x m )e lwei x \ 
e iu,0 2 .x mj 6» 2 ■T(x m )e l " e2 ' x , 6> 2 -r 7 (x m )e l " e2 - x 

V e iu0 f*™, 9 Q ■ r(x ro )e^°«- x , Q • r7(x m )e^ e « x / 



The multi-frequency imaging algorithm developed in 0. EH EH Il9ll is given hereafter; for several fre- 
quencies {u>k '■ k = 1, 2, • • ■ , 7^}, Singular Value Decomposition (SVD) of A(lu) is performed as follows 



ak)=uk)§K)v*K) » f™(wOUfflK)v;(w k ) 



where (J m (u!k) are nonzero singular values, and U m (u;fc) and V m (ojfc) are respectively the left and right 
singular vectors of A(wfc). Then, for a suitable vector a e K 3 \{0} and z e Q, define vectors E(z; u^) € C Px 1 
andF(z;w) e C Qxl as 



E(z; Wfc ) = ( a ■ (1, flOe"^*, a • (1, tf 2 )e"^' z , • • • , a • (1, tf P )< 



— IWtfp -Z 



F(z; Wfc ) = (a • (1, ^Je*"*-, a • (1, 9 2 )e 1 ^, • • • , a • (1, Q )e^ e « z 
respectively. Consequently the corresponding unit vector can be generated: 
Wj(z;w fc ):=— — <- and W L (z;w fe ):= 



|E(z; Wfe )l ' ' ' |F(z;w fc )| 

Subsequently, we consider the following multi-frequency imaging function (see HQ EE H) 



K M k 



k—1 m—1 



J2 W}(z; Wfc )U m K) Wl(z; Wt )V m K) 



where a denotes the complex conjugate of a e C. Then, by virtue in JHEHEH-EHl/ we get 

Z(z,tf)~x(T), 



(4) 



(5) 



where A ^ B means that there exists a constant C such that A = CB. 

Although (O indicates why the proposed imaging function produces an image of thin inclusion(s), a 
mathematical analysis is still needed for explaining phenomena such as the appearance of unexpected 
ghost replicas. In the following section, we analyze the imaging function 10 and improve it for better 
imaging performance. 



3. Introduction to weighted multi-frequency imaging function 

3.1. Analysis of multi-frequency imaging function (fSJ) 

Now, we analyze the multi-frequency imaging function fl5). To do this, we assume that P is sufficiently 
large and number of nonzero singular values Mk is almost equal to M, for k = 1, 2, • • • ,K. We then obtain 
the following results. 
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Theorem 3.1. Assuming that K is sufficiently large, Q becomes 

M Q 

EE 



i 



, Q is small and luk — > +00, 



X(z,K) ~ < 



1 = 1 g=l yj | Z - X m | 2 - (0 q ■ (Z - X m ))' 

M 



u) K - wi 



"777" E ^*(l z - x mh w K) - *(l z - x m|; w i)j 7 Q is large and uj k < +00, 
S(z, x m ) , Q zs Zarge and wj{ — > +00. 



Here, the function w) zs graer; fey 

:= Lo(j (Ljt) 2 + Ji(o;i) 2 ^ +^ Ji(xt) 2 dx, 

where J v (t) denotes the Bessel function of order v and of the first kind and 6 denotes the Dirac delta function. 

Proof. First, assume that Q is sufficiently large and that lu k < +00. We denote Ai9 p := \d p - 
p = 2, 3, ■ • • , P, and #1 := — i9jv|- Notation A9 q can be defined analogously. 
Since the following property (see [4]) holds 

W, 7 (z; ui k ) ~ U m (w fc ) and W L (z; w fc ) ~ V m (w fe ), 

plugging the above relationship into © and applying fiol Lemma 4.1], we observe that 

K M 



X(z,K) 



E E (w}(z; Wfc )W, 7 (x m ; Wfc )") [w2(z; Wfc )W L (x m ; Wfc )) 
fc=l m=l ^ ' 
K M / P 

EE(E^ p - (z - Xm) 



fc=l m=l v p=l 



x 2tt 

fc— 1 m— 1 v p— 1 
if M 



1 

4^2 



1\ iVJ y „ 

E E ( y e^-^-^de 

Then, applying the indefinite integral (see ||23L Page 35^) 



8=1 



9 =1 
1 



E E (^^e^'V^-^Atfp) ( _L^ e -^ 9 .(z-x„ l)A0 

^2 / Jo(^|z - x m |) 2 rfw 



JoW^dt = t J (zT + Ji(t)^ + / Jiftrdt 



with a change of variable cj|z — x m | = t yields 

M 



X(z,K) 



E 

in— 1 
M 

E 



1 



UJK Z — X T) 



z x m| Jwi|z-x m | 



J (t) 2 dt 



\) 2 duj 



uj k - uii 



Jo(u K \z - x m |) 2 + Ji(ujk\z - x m |) 2 + / Ji(w|z-x r , 

Jo 

J (wi|z - x m |) 2 + Ji(wi|z - x m |) 2 ^ + J Jx{ui\z - x m \) 2 dwj 



^ ($(|z-x m |;o;A')-$(|z-x m |;a;i) 



1 Iti our knowledge, there is no finite representation of the integral J ,J v {t) 2 dt for an integer u. 



Consequently, if we assume ljk — > +00, then 

I(z,K) ~ 5(z,x m ). 

Next, we assume that Q is small and that cjk — > +00. Then 

K M / P \ / Q 



l(z,K) 



m— 1 ^ p— 1 

J2 / e 4 ^- (z - Xm) JoMz - x m |)d, 
m= i „=i 1 



iw fc e,-(z-x TO ) 



Since 0, e S 1 , we note that 



|z - x m | 2 - (0 q ■ (z - x m )) = |z-x,, 



1- 



> 0. 



Hence, applying the identity (see |9, formula 6.671, page 717]) 



yields 



r°° 1 
/ e lat JJbt)dt = - cos ( sin" 1 - 

Jo Vb 2 - a 2 [ V b 

e ^.(z-x ro ) Jo(w | z _ Xm | )2dw 



. / . _i a 
3in ^ sm — 
V b 



for a < 6 



lim 



Z - X m | 2 - (flg • (Z - X m ))" 



□ 



Remark 3.2. Theorem 1331 indicates certain properties 0/©, u?/zic/z are summarized as follows: 



1. Jo(x) ftas an oscillation property. Hence, if a small value of K is considered, the identification of shape o/T 
zt>z7Z fee interrupted by ghost replicas. 

2. On the other hand, if a large number of frequencies are applied, X(z, K) will exhibit an accurate shape o/T. 
Note that this fact can be validated via Statistical Hypothesis Testing (refer to [4] for more details). 

3. When Q is small, J(z, K) plots a large magnitude at z satisfying 

z = x m e T and q = ±- —. 



This means that I(z, K) produces not only the shape o/T but also unexpected ghost replicas. 
4. When low frequencies ujk are applied such that ujk\z — x m | w 0, 1(z, K) will fail to produce the shape o/T. 
However, if sufficiently high frequencies ujk are applied, the term 



0(|z - x m |,w) := / Ji(w|z - x m |) 2 du; < 



0(lo k ) 



becomes negligible (see and hence the shape o/T can be easily constructed by l(z, K). This is the reason 
why applying high frequency yields good results. 
5. For a small value of K, ifujK — > +00 then it is expected thatl(z, K) will yields a good result, however, this 
is an ideal assumption. 
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3.2. Permeability contrast case: why do two ghost replicas appear in the numerical experiments? 

In section |2j although we did not focused on the vector a, it is important. In order to determine the 
influence of a on imaging performance, let us define vectors E(z; Wfc) € C Pxl and F(z; ui) € C^ xl in ((Hi as 



E(z; w fc ) = ^ e - iu,,?I - z , e"*"**-, . . . , e -i"*p-*j } F (z; Wjk ) = ^j*"® 1- ", e iw ° 2 ' z , • • • , e lwe « z J , (7) 
and corresponding unit vectors Wj(z;ujk) and Wi(z;u;fe). Then applying jiol Lemma 4.1], X(z,K) be- 



comes 



X(z,A-) 



k=l m=l ^ ' ^ ' 

K M / P \ / Q ^ 

E £ £ *P • ( x ( x »0 + rKx™))^^ 2 -^ £ V (r(x m ) + J7 (x m )) e ^ e '-( z - x '"' 

fe=l m — 1 \p— 1 / \q—l / 

K M / p \ 2 

E E ( / ' ( T ( x ™) + r?(x m ))e^ e '( z - x -)d0 ) 

/"(Lux / \ 2 ^ /"Wft: 

^ / f (z - x m ) ■ (r(x m ) + T7(x m ))Ji(w|z - x m |) j duj ~ ^ / Ji(w|z - x m |) 2 dw. 



Note that J\(u)x) is maximum at two points xi and £2, and is symmetric with respect to x (see FigureQJ. 
This means that l(z, K) plots on T and two ghost replicas in the neighborhood of T (see Figure|2]l. This 
means a must be a linear combination of r(x m ) and 77 (x m ). Unfortunately, we have no prior information 
of T. Hence, finding an optimal a. is an interesting research topic. 
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Figure 1: Graph of y = Ji(u>\x\) 2 for ui = 27r/0.5. 



3.3. Weighted multi-frequency imaging function: Introduction and analysis 

Remark l3.2l shows the optimal conditions (P and Q are sufficiently large and luk — > +00) for imaging. 
However, real-world applications, a high frequency such as +00 cannot be applied. Therefore, in order to 
improve 10, we propose a multi-frequency imaging function weighted by (u!k) n as follows: 



I(z, K; n) 



K M k 

EE' 

fe=l m=l 



(a; fc )™ W}(z;w fc )U m (w fc ) Wl(z; Wfc )V m ( Wfe 



(8) 



Note thatl(z, K; 0) = I(z, if). Consequently, the following question arises: 
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Does a large value of order n contributes significantly to the imaging performance? 
The answer is no. To find out how, the structure of © is explored as follows. 

Theorem 3.3. Let K and wk have sufficiently large values and lok < +00. Then for a natural number n, © 
satisfies 

I(z,K;ri)~ (*(l z - x m|,WK;n)-$(|z-x ro |,u;i;ra)), (9) 



where $ is given by 



<$>{t,u-n) J (ujty + Ji(wt) 



Proof. First, let us assume that n is an odd number, say n = 2 s + 1. In order to examine the structure of (8), 
we recall a recurrence formula (see I23I page 14]) 



t 2n+i j Q (tydt = 



1 



2n + 1 



2n+2 



Jo{ty + Ji(t) 2 )+n H 2n J (ty+nt 2n+1 J (t)J 1 (t)~2n- i / t 2n - x J (t) 2 dt 



Then, applying dTOb with a change of variable 10 \ z — x m | = t in (|8]l induces 



(10) 



X(z, if; 2s + 1) 



M 



z x |2s+2 



t 2s+L j {tydt 



Y I *i(|z-x m |,WK,cJi;s) + * 2 (|z-x m |,WK,wi;s) 



* 3 (|z - x m |,WK,wi; s) - * 4 (|z - x^ljWjfjWijs' 



where 



*i(t,wjr,a;i; s) 
i &2(t,UK,u 1 ; s) 



2s+2 



(2s + l)i 2 
s(^) 2s+1 



4s + 2 

2s 7 / . .\2 



Jo(wit) 2 + Ji(wit) : 



and 



(2s + l)t 
* 4 (|z - x m |,Wif,tJi; s) := 



s(cji) 2s+1 

Jo(uJ K t)Ji(uJKt) - — — — Jo{wit)Ji(uit), 



(11) 



(2s + l)t 
2s 3 



(2s + l)|z-x m | 2 



l{z,K;2s-l). 



Now, we assume that ojk is sufficiently large, z satisfies z 7^ x m and |z — x m | ^> 4^7—- Then, applying 



the asymptotic form of the Bessel function 



J v {t) ~ \ — cos t - 

7TE 



yields 



CJff Z - XmpT 



COS Wif z — x r , 



Z/7T 7T 

T ~ 4 



4 << ^7^ C0S K |z_Xml ~i;' 



(12) 



and 



J (uj k \2 - x m |) J\{u K \z - x m |) 



UJ K \Z - X m |^7T 



cos I cj k \z - x m \ — -j I cos I ojk\ z — x m | - 



3?r 



< Q COS I Wif I Z - X TO | —I COS I Wjf|z - X m | - — 

Hence, we can observe that 

* 2 (|z -x m |,w K ,a;i;s) < O ((u k ) 2s+2 ) = *i(|z - x^wj^wi,; s), 
^(Iz-Xrol^K.WlJs) < ((w fc ) 2s+2 ) = *i(|z -x m |,cjA-,Wi;s), 

* 4 (|z - x m |,a,K-,wi; s) « ^^^ (z, K; 2s - 1) = O (K) 2s+1 ) . 

This means that the terms ^2, ^3, and ^4 are dominated by $1. Hence, we conclude that for an odd 
number n, © holds. 

Next, let us assume that n is an even number, say n = 2s. Similar to the case of an odd number n, we 
recall a recurrence formula (see [23, page 35]) 



Jt 2n J Q (t) 2 dt = A n {t)J (t) 2 + B n {t)Mt),h(t) +C n {t)J 1 {t) 2 + D„ J Mtfdt, 



(13) 



where 



n+l 



A n (t) = drt 2r -\ B n {t) = J2 brt 2r , C„(t) = £ c n t 2r +\ 



r=l 



and a r ,b r ,c r , and £>„ are constants. Then, applying a change of variable t = ui\z — x m | in J131 yields 

M , /-WKl»-Xm| 



X(z,if;2s) 



1 



where 



m— 1 

I * 5 (|z - x m |, wk, wx; s ) + *e(|z - x TO |,WAr,a;i;s) 



+ *7(|z - x m \,u K ,LJi; s) - * 8 (|z - x m |,Wif ,wi; s) 

1 s+1 / \ 
^(i^A-^iss) := (^ar(wA-t) 2r " 1 J (w A 't) 2 - a r (wii) 2r ~ 1 J (a;it) 2 J, 

(b r (u K t) 2r jQ(u K t)Ji(LU K t) -b r (uj 1 t) 2r Jo(uj 1 t)J 1 (uj 1 t)j, 



*7(i,W}f,wi; s 



and 



t 2s+l 



* 8 (|z - x„,|,us:,wi; s) 



D n 



-l(z,K). 



Now, we assume that u)k is sufficiently large, z satisfies z 7^ x m , and |z — x m | 3> 4^—- Then applying 
(O yields 

2^a r (w A '|z - x m |) ^Z^M-H <max s {a r } — ^ . 

r— 1 r— 1 ^ ' 
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Accordingly, we obtain 

4- 5 (|z-x ro |,u> K ,wi;s) ~ [ (w^) 2s+1 Jo(w/f|z-x m |) 2 - (wO^+VqKIz-x™!) 2 



* 6 (|z-x m |,o;K,wi;s) < ^ 5 (|z - x m |,w_ft:,a;i; s), 

* 7 (|z-x ro |,w K ,Wi;s) - ^(w K ) 2s+1 Ji(w /f |z-x m |) 2 - (u^+ViMz-x™!) 2 

* 8 (|z-x m |,^,^i;s)« D n l(z,K) -0(K) 2i+1 )- 

Hence, we conclude that for an even number n, |9]l holds. □ 
Remark 3.4. Based on Theorem \33[ we can determine following properties of®, as summarized below: 

(PI). Increasing n in © does not influence the imaging performance. In fact, n = 1 is the best choice for obtain- 
ing good results because the terms *;(|z — x m |, ujk , s), I = 2,3, • • • , 8, in ® are eliminated. 

(P2). Similar to Q, applying low frequencies ujk such that u>k\z — x m | rs degrades the results. Hence, sufficiently 
high frequencies must be applied. 

(P3). When n is an odd number, l(z, K; n) does not contain 



/>OJK 

/ Ji(w|z - x m \) 2 du 



but it is a part o/I(z, K; n) when n is an even number. This means that, although l(z, K; n) offers good 
results, large numbers of unexpected replicas will appear in the map ofl(z, K; n) when n is even and ujk is 
large. Hence, for an integer s, Z(z, K; 2s — 1) offers better results than Z(z, K; 2s). 

4. Numerical experiments 

In this section, a set of numerical experiments are performed. Throughout this section, the homoge- 
neous domain £1 is selected as a unit circle and two jj are chosen to describe thin inclusions T j as 

7! = {{z - 0.2, -0.5z 2 + 0.4) : z G [-0.5, 0.5]} 
72 = {(z + 0.2, z 3 + z 2 - 0.5) : z € [-0.5, 0.5]}. 

The thickness of all the thin inclusion(s) is set to h = 0.015, and £q = Mo = 1 an d £j = (J-j = 5. For given 
wavelengths \k, the applied frequencies are ujk = f^v which are varied between ui\ = ^ and ujk — §3 • In 
every example, K = 10 different frequencies are applied. 

A set of vectors : j = 1, 2, • • • , P} and {6i : I = 1, 2, • • • , Q} on S 1 are selected as 

, 2nj . 2nj\ , „ / 2ttZ . 2ttZ 
Wj = I cos -^r- , sin I and &i = — | cos , sm 



P P ) \ Q Q 

respectively. In our examples, we set, P = 24 and Q = 20. Vector a in |(4) is set to a. = (1, 0, 0). 

For every example, white Gaussian noise with a Signal-to-Noise Ratio (SNR) of lOdB is added to the 
unperturbed boundary measurement data using the MATLAB command awgn in order to demonstrate the 
effectiveness of the proposed algorithm. For discriminating non-zero singular values, a 0.01— threshold 
strategy is adopted (see [19] for instance). 

First, let us consider the imaging of Ti when only magnetic permeability contrast exists, at a fixed fre- 
quency w = Itj. From the result in section l3l2l it is clear that two ghost replicas appear in the neighborhood 
of Ti when we select vectors E(z; Wk) and F(z; w) in (0. Note that based on the shape of Ti, r(x m ) « 
and rj(x m ) ps 1. Hence setting of a = (0, 1) yields a good result while a = (1, 0) offers a poor result, as 
shown in Figure 12 
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Figure 2: [Permeability contrast case only] Maps of I(z, 0) with vectors in (7) (left), a = (1, 0) (center), and a = (0, 1) (right) when 
the thin inclusion is Ti . 
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Figure 3: [Permeability contrast case only] Maps of Z(z, n) with a = (0, 1) for n = 3 (left), n = 5 (center), and n = 10 (right) when 
the thin inclusion is Ti . 



Figure [3] shows the effect of the number of frequencies on the imaging performance. By comparing the 
results with those in Figured we observe that applying multiple frequencies guarantees better results than 
applying a single frequency does. 

Hereafter, we consider the imaging of thin inclusions when both permittivity and permeability contrast 
exist. Figure S] shows the maps of I(z, 10;n) for n = 0,1, and 2. The overall results show that although 
I(z, 10; n) offers good imaging performance, odd number of n yields better results than an even number of 
n does, refer to |(P3)| of Remark l3~4l 
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Figure 4: Maps of I(z, 10; n) for n = (left), n = 1 (center), and n = 2 (right) when the thin inclusion is Ti . 



Now, we consider the maps of l(z, 10; n) for n = 5,6, and 7. The corresponding results illustrated in 
Figure [5] indicate that increasing n not only contributes negligibly to the imaging performance but also 
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generates large numbers of unexpected replicas. This result supports (PI) of Remark l3.4l 
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Figure 5: Maps of X(z, 10; n) for n = 5 (left), n = 6 (center), and n = 7 (right) when the thin inclusion is Ti . 

Figure [6] shows the maps of l(z, 10; n) for n = 0,1, and 2 when the thin inclusion is T2. Similar to 
Figure|H the map of I(z, 10; 1) demonstrates very good imaging performance but the resolution is poor at 
z = (0.6 — 0.24) on the large curvature. 
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Figure 6: Same as Figure|4]except that the thin inclusion is T2 . 



For the final example, we consider imaging of multiple inclusions Ti U T2 with e-y = /zi = 5 and 
£2 = fJ-2 = 10. Figure [7] shows the maps of l(z, 10; n) for n = 0,1, and 2. Similar to the previous examples, 
we can see that the map of I(z, 10; 1) yields the most accurate shape of Ti U T 2 , but because Ti has a much 
smaller values of permittivity and permeability than T2, 1(z, 10; n) plots Ti much smaller magnitude than 
T 2 . 
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Figure 7: Same as Figure|4]except that the thin inclusions are Ti UT2. 
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5. Conclusion and perspectives 

We introduced a multi-frequency algorithm weighted by applied frequency for imaging thin electro- 
magnetic inclusions and discussed its properties. The structure of the imaging function and the corre- 
sponding numerical experiments show that the proposed imaging algorithm is effective and robust against 
a large amount of random noise; however application of a large number n in ((8) does not guarantee a good 
imaging result. 

We considered a full- view inverse scattering problem, but based on our contributions jl^,ll8tl , the pro- 
posed algorithm is also applicable to the limited-view inverse problem when the measured boundary data 
is affected by randomly distributed small inclusions. Therefore, extension of the algorithm to the limited- 
view inverse problem affected by random scatterers will be an interesting research topic. 
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